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ABSTRACT 

Wc hypothesise that planets are made by tidal downsizing of migrating giant planet 
embryos. The proposed scheme for planet formation consists of these steps: (i) a mas- 
sive young protoplanetary disc fragments at i? ~ several tens to hundreds of AU on 
gaseous clumps with masses of a few Jupiter masses; (ii) the clumps cool and contract, 
and simulta neously migrate closer in to the parent star; (iii) as earlier suggested by 
iBossI ()1998l ). dust sediments inside the gas clumps to form terrestrial mass solid cores; 
(iv) if the solid core becomes more massive than ^ 10 M^, a massive gas atmosphere 
collapses onto the solid core; (v) when the gas clumps reach the inner few AU from 
the star, tidal shear and evaporation due to stellar irradiation peel off the outer metal- 
poor envelope of the clump. If tidal disruption occurs quickly, while the system is still 
in stage (iii), a terrestrial planet core is left. If it happens later, in stage (iv), a metal 
rich gas giant planet with a solid core emerges from the envelope. 



1 INTRODUCTION 



There are currently two competing theories for planet for- 
mation in which the planet making action starts fr om oppo- 
site en d s. The core ac c retion ( CA) model (e.g., see lSafronovl 



(|1969D : iPoUacket af] (1 19961 ): llda fc LinI (|2008l ). and chap- 
ters 4-6 in lArmitagel (|2010l )) stipulates that planets form in 
a bottom-up scenario. First, microscopic dust in the pro- 
toplanetary disc grows into grains. Th ere is then a po orly 



understood but necessary step (e.g., Wetherilll 199d) of 



grains joining up som ehow (see lYoudin fc GoodmanI uOOSi : 



Ijohansen et al.l 120071 ) to make km-sized objects (planetesi- 
mals) , which then collide and stick to form ever larger solids 
up to terrestrial planet masses. When the solid core reaches 
a critical mass of about 10 Earth masses, a massive gas en- 
velope builds up by accretion of gas from the disc. 

In contrast, the Gravitational Instability (GI) model 
suggests that giant planets form "big", e.g., directly by con- 
traction of gaseous c ondensations born in a massive self- 
gravitating disc (e.g.. iBodenheimeij [TOTJ : IBossII1997| ). Gi- 
ant planets in the Solar Syste m contain massive solid cores 
IjFortnev fc Nettelmannll2009l ). as pr edicted by t he core ac- 
cretion model. In the context of GI, IBossI ([1993) suggested 
that such cores could r esult from dust sed imentation inside 
gaseous proto-planets. IBoss et al.l (|2002l ) showed that re- 
moval of the gas envelope by photoionisation due to a nearby 
OB star could make the ice giants of the Solar System. 

Of these two models, the GI is the one that received 
the lion's share of the criticism. The most notable prob- 
lems of the model are: (a) protoplanetary discs cannot 
form giant pl anet embryos at the location of Jupiter (e.g.. 



piter (e 

Rafikovll2005h cont r ary to e arlier results b y IBossI (j 19971 ): (b) 



Helled fc Schubert! (|2008l ): iHelled et all (|2008l ) found that 



dust sedimentation is too slow a process to yield observed 
solid cores in giant embryos of mass ^ Jupiter's mass; (c) 
OB stars are too rare to explain the abundance of Jupiter- 
like planets observed in extrasolar planetary systems, and 
(d) there seems to be no way to form terrestrial planets like 
the Earth. 

The goal of our Letter is to show that this criticism 
has long outlived itself. All of the points raised above can 
be addressed if one upgrades the 1990-ies version of the 
GI model with more modern ideas - formation of Giant 
planet Embryos (GE hereaft er) far from t he parent star 
where the process is allowed (|Rafikovll2005l ). and then let- 
ting them migrate radially i nwards, as massive planets do 
(|Goldreich fc Tremain"3ll980l ). 

We have recently revisited (paper I, Nay akshin 2010a, 
to appear in M NRAS) the suggestions of IBossI (|l998t ): 
IBoss et al.l (|2002l ) on dust sedimentation in GEs assuming 
that they are formed at Rp ~ 100 AU. Such GEs are cooler 
by an order of magnitude than those studied by Helled 
and co-authors. GEs with masses below about 10-20 Jupiter 
masses (Mj) are found (paper II, Nayakshin 2010b, submit- 
ted to MNRAS) to be excellent sites of solid core formation. 
Thus addressing point (a) addresses point (b) automatically. 

Below we estimate the migration time of the embryos 
in a realistic young massive protoplanetary disc. We note 
that as embryos migrate closer to the parent star, their Hill 
radius decreases faster than they can shrink due to internal 
cooling. When these two are equal, we assume that the GE 
can be disrupted by the tidal field of the star, which nor- 
mally occurs at a few AU distance. However, we find that 
grain sedimentation time is shorter than the migration time 
for plausible parameter values. Thus, by the time the em- 
bryos are due to be disrupted, they may already contain 
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massive solid cores. Depending on the properties of the GE 
at disruption, removal of the outer metal poor gas shell re- 
veals either a terrestrial planet or a metal-rich giant planet. 
We therefore argue that such a "tidal downsizing" scheme 
may expl ain planets of a l l types without recourse to plan- 
etesimals (|Safronovlll9"69l : IWetheriljll990l ) . 

Our hypothesis is a hybrid scenario: giant planet forma- 
tion starts with dust s edimentation inside a self-gravitating 
gas clump (jBossll 19981 ) , as in the GI model, but continues as 
in the CA model with accretion of a massive metal-rich gas 
envelope from within the GE. Additionally, solid cores left 
by the GE disruption in the protoplanetary disc may grow 
further like in the CA model - by accreting smaller solid 
debris from the disc. Giant planets may accrete more gas 
and also migrate radially, as in the "standard theory" . 



2 GIANT PLANET EMBRYO MODEL 

2.1 Birth and contraction 

Massive young gas discs are gravitationally unstable and 
fragment on clumps when the disc cooling time is no longer 
than a few t imes the local dynamical time, l/Jl (e.g., 
iGammid [JOoH '). where O. = {GM^/Riy^^ is the Keple- 
rian angular frequency at radius Rp from the star of mass 
M*. This l imits the fragmentation region to Rj, ^ 30 — 100 
AU fe.g . iRafikovl l2005l : IStamateUos fc WhitworthI l2008l : 
iMeru fc BatdbOld ). At the point of marginal gravitational 



stability the density in the disc midplane is about the tidal 
density at that point, pt = M*/(27r_Rp). Thus, when the disc 
just fragments, the disc cooling time is about the free-fall 
timescale for the material in the disc, iff ~ l/\/Gp. As the 
clumps contract further, their cooling time becomes longer 
than their free-fall times. T he initial state of a G E is thus 
that of a polytropic clump (|Cameron et al.lll982l '). similar 
to the "first cores" studied in the context o f star formation 
IjLarsonl Il969l : iMasunaga fc Inutsukal [200y) . For simplicity 
we assume that GEs evolve at constant mass after their for- 
mation. 

Introducing dust opacity as a power-law in temperature 
and assuming that the embryo is optically thick, we obtained 
an approximate analytic solution for the radiative cooling 
of the giant embryo in paper II. For example, for k{T) = 
K.o{T/lOK), we have 



T (t) = To 



1 + 2 



1/2 



146 mt/\,-*/' 



1 + 2 



1/2 



(1) 



where To is the temperature of the embryo at formation, 
with mi = Mcmb/lOMj being the mass of the gas clump in 
units of lOMj, kt — kq/O.OI, and to is the initial cooling 
time of the clump: 



, oon 2/3 1/9 

to — 380 years ttIj kJ 



(2) 



Note that the cooling time of the embryo increases with 
time, so that tcooi(i) oc t when t ^ to. At a, constant mass 
the embryo radius decreases with time as 



Jtcinb — 



-Ro 



[1 + 2i/io 



1/2 



= 6.0 AU 



4/9 



'^'To 



-1/2 



(3) 



At large t, independently of Memb, all GEs contract as 



Rr. 



(t) = 0.8 AU K^^ 



lO'' yrs 



1/2 



(4) 



2.2 Dust growth and sedimentation 

We echo the approach of IBossI ([1993) here; the resulting 
grain growth time scale in a constant density embryo, for a 
grain with initial size ao to reach the final size a is 

igi-.o = f „%n In — , (5) 

-KjgpoGRo ao 

where po = 3Momb/47r_Rg^bi Ca and fg ^ 0.01 are the initial 
embryo mean density, the sound speed, and the grain frac- 
tion by mass, respectively. In paper I we improved on this 
model by allowing the embryo to contract, in which case the 
growth time is a geometric combination of the cooling and 
the initial grain growth time, igr = tg tpi o- 



^f^'fzTkV'' ( 



ln(a/ao) 
20 



4/7 



(6) 



where /_2 = /g/0.01. The vaporisation time, tvap, is the 
time it takes for the embryo to heat up to the vaporisation 
temperature of Tvap ~ 1400 K. This is trivially obtained 
by solving the equation T(tvap) = Tvap (using eq. [l|, and a 
good approximation is 



= 1.5 X 10 years 



T 



1400 K 



(7) 



Grains sediment to the centre of the core when tgr < tvap, 
and vaporise before they could sediment if fgr > tvap. 

The mass of refractory material (e.g., silicates) that can 
be used to build a massive core is expected to be about ~ 1/3 
of the total heavy el ements mass, or about 20 Earth masses 
for Memb = lOMj (|Boss et all 120021 '). However, due to en- 
ergy release associated with the solid core accretion, strong 
convective motions near the core may develop. The grains 
also melt at too high core accretion rates. Thus the final 
mass of the core is usually much smaller than the maximum 
condensible mass (paper II). 



2.3 Radial migration of embryos 

Gravitational torques betw een the disc and massive 
planets are signifi cant (e.g., iGoldreich fc Tremaind Il980l : 
iTanakaet alll2002l '). While GEs are smaller than their Hill's 
radii (see below), we assume that they migrate at similar 
rates. Under this assumption, embryos migrate via type I 
regime if Mcmb ^ Mt = 2Mt{H/R)'^ and via typ e II regime 
if Me mb > Mt, where Mt is the transition mass (jBate et al.l 
[2003) • For a self-gravitating, Q — 1 disc, the type II migra- 
tion time is 



,(8) 



J- i-r, \ -lo-l ^P A in4 -1 o3/2 /0.2itp \ 

tu{Rp)^a n — 1=4x10 yrsQ_i7?2 ( — Tf^ ) 

where a_i — a/0.1 is the disc viscosity parameter and R2 — 
i?p/100AU. Type I migration is faster: 



ti (Rp) 



D.- 



M« H 

Memb Rp 



3 X 10"^ yrs R. 



3/2 



H 



0.2Rp 



(9) 



We can improve these estimates by building a simple disc 
model that would specify how disc properties change with 
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R. We define the stellar mass doubling time scale as t^b ~ 
M^/M where M is the mass accretion rate on the star, as- 
sumed to be constant in time and radius inside the disc. By 
the order of magnitude, f db should be comparable to the free 
fall time of the host gaseous cloud from which the star forms, 
e.g., ~ 10^ yrs for a 1 Mp ) cloud at th e typical interstellar 
temperature of 10 K (e.g.. LarsonI igGg). 

Using the standard IShakura fc SunvaevI (|l973l ) accre- 
tion disc formalism, we first find H{R) and Md{R) in the 
innermo s t disc . We then find the radius beyond which the 
iToomrd (| 19641 ) parameter falls below unity. In the outer 
sel f-gravitatin g regio n of the disc we follow the treatment 
of iGoodmanI ((2003) for a disc with no internal energy 
sources and enforce Q ^ {H/R) {M,/Md) = 1. Finally, 
for a constant accretion rate disc, M — SivaH^fl'Ed, where 
Ed « Md{R)/{'KR?) is the column depth of the disc, we can 
solve for 



H_ 



Safiid 



1/3 



(10) 



The accretion rate disc model built in this way is schematic 
but correctly predicts the expected location of the self- 
gravitational region at Rp > few tens of AU to 100 AU. 
Using this disc model we are now able to calculate the mi- 
gration rate of a GE at an arbitrary R^ for a given tg^. We 
follow the fit to numerical results by iBate et al.l (|2003l ) to 
smoothly join the type I and type II migration regimes. The 
fastest migration then occurs for embryo mass Afcmb = Mt- 

2.4 Tidal disruption 

The Hill's radius of the embryo is (for Af* = IM©) 



-Rh ^Rp 



Me, 



3M, 



1/3 



= 0.15 RpTn^ 



1/3 



(11) 



In general, radial migration accelerates as the GE moves in, 
whereas the contraction of the envelope slows down with 
time (cf. equation|4]). Therefore, as planets migrate inwards, 
there is a point where the embryo's radius 7?emb becomes 
comparable to the GE Hill's radius. The tidal field from the 
parent star at this point starts to peel off the outer layers 
of the embryo. The tidal radius, Rt, is defined as the radius 
where i?cmb = fltRn {rjt ^ 1). Using the type II migration 
estimate as an example, we insert eq. [8]into eq. |4l and find 

H 



Rt = 2.8 AU %- alYi?, 



l/2p-3/4 1/2 



'1/3 



0.2R 



(12) 



In Appendix we show that irradiation by the parent star can 
heat up the outer layers of the envelope and even disrupt 
the whole GE. The effect appears less important than tidal 
disruption for Solar type stars, but may become dominant 
for higher mass stars. 



3 AN ILLUSTRATIVE EXAMPLE 

Figure [T] shows the results of our giant embryo model de- 
scribed in §2 for embryos of 3, 6 and 10 Jupiter masses, 
shown with different colours. The parameters used for the 
figure are tdb = 10^ yrs, a = 0.1, kq = 0.01, M* = 0.5 M© 
(as the star is assumed to be in the midst of its growth, 
rather than close to being completely assembled) , Solar met- 
alicity, and the starting radial position is Rp = 100 AU. 
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Figure 1. Evohition of GEs of mass 3, 6 and 10 Jupiter masses 
(black solid, red dotted and blue dashed, respectively) calculated 
as described in §2. Upper panel: Radial position of the embryos 
versus time. The GEs start off at /J = 100 AU. The coloured 



asterisks and text "core formation" show t. 



the time when a 



solid core should form inside the embryo, with the colour matched 
to that of the curve. Midcile: Radius of the first core (power-law 
like curves), and the tidal radius at the current location of the 
first curve. The first core is assumed to be completely disrupted 
when the two sets of curves meet, i.e., when the tidal radius drops 
below -Rcmb- Lower: Temperature of the first cores as a function 
of time, with core formation times marked as above. 



The upper panel shows radial migration of the embryos. 
Note that initially all migrate at about the same rate (in the 
type II regime) , but then the least massive embryo overtakes 
the other two. The least massive one (black solid curves) 
migrates the fastest as it happens to fall in the minimum of 
the migration time curve, betwe en type I and type II regimes 
(see fig. 11 in IBate et al.l 120031 ). The asterisks and the text 
"core formation" mark the time t = igr, when the solid core 
is assumed to form in the centre of the GEs (cf. eq. |S|. 

The middle panel shows the embryo size, i?omb, and 
the respective Hill's radius, for the three cases considered. 
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We assume that when the two sets of curves intersect, the 
embryo is disrupted, leaving only the core, which migrates 
at a negligible rate during the calculation. This is why the 
radial position curves, Rp{t), approach a constant value in 
the upper panel at late times. 

The lower panel shows the respective temperature evo- 
lution of the embryos (cf. eq. [T]). The least massive embryo 
forms the core last, as it is the least dense of the three. How- 
ever, it is also the coolest one. Due to this, the least massive 
embryos may lock most of their "metal" content into grains 
and actually yield heavier solid cores (see paper II). On the 
other hand, the least massive embryo is disrupted sooner (in 
accord with ea. ll2p . at i?p ~ 6 AU. 



4 THE TIDAL DOWNSIZING HYPOTHESIS 

4.1 Planets: leftovers of disrupted GE ? 

We found that for plausible conditions, the GE initial cool- 
ing time, the dust growth time, the grain vaporisation time 
and the GE migration time are related by 



< tn 



(13) 



respectively. This order of timescales implies that (1) the 
GE contracts quickly initially, which protects it from an im- 
mediate tidal disruption at the point where it is born; (2) 
grains have enough time to grow and sediment to the centre 
of the GE, forming a massive solid core there; (3) since the 
cooling time becomes asymptotically long f ^2.1[) . the embryo 
is tidally disrupted at the inner few AU. 



4-. 1.1 Terrestrial planets 

If the migration time is shorter than the vaporisation time, 



t„ 



< iv 



(14) 



and the solid core mass, Mcorc, is smaller than the critical 
core mass. Ma ~ 10 M® (see below), the gaseous component 
of the GE is almost completely disrupted, e.g., disassociated 
from the solid core. As the density of the solid cores is a few 
g cm~^, the solid cores and any smaller debris (if present) 
are not disrupted, and instead are released into the proto- 
planetary disc around the GE location. The least massive of 
the debris bodies may remain at that location in the disc for 
as long as xi lO" years, when the gaseous disc is presumably 
dissipated. After the disruption, low mass solid cores are en- 
dowed with a tenuous gas atmosphere (paper II), much of 
which may be lost thereafter. 

4- 1-2 Giant planets with solid cores 

In paper II we argued that when the solid core grows more 
massive than a critical core mass, Ma, the atmosphere near 
the core becomes so dense that it is gravitationally unstable 
and starts to collapse on the solid core. This "core-assisted" 
collapse is quite similar to that found in the classical CA 
model fe.g.. lMizunolll98d : [Follack et al.lll996l 'l. except that 
gas is accreted on the core from the GE rather than the disc. 
A proper numerical calculation is required to find the exact 
value of the critical mass, Ma, but one can limit it from be- 
low by assuming an isothermal gaseous envelope (e.g., §6.1 



in lArmitageil2010l ). which yields Afcr ~ 3M0 (see paper II). 
The maximum value of Afcr can be obtained by employ- 
ing the classical CA model results - the r adiative zero so- 
lution (|Stevensonlll982l : llkoma et al.ll2000l ). and this yields 

Mcr ^20Mg). 

We would thus posit that, if the tidal disruption of the 
GE occurs when the core mass exceeds Afcr, the outer metal- 
poor low density layers of the GE are peeled off, but the 
inner high density solid core plus the bound gas atmosphere 
remains. The disruption leftover must be a metal rich (com- 
pared to the parent star) giant planet, as are the giant plan- 
ets in the Solar System. 



4.1.3 Giant planets without solid cores? 

The GE may go through the second wave of collapse when 
the temperature reaches about 2000-2500 K (|LarsorJll969l : 
lBodenheimenll974l : iMasunaga fc Inutsukall200Cl ) which can 
compact the gas component to densities as high as ~ g cm~'^. 
The time to reach that temperature is at least a few times 
the vaporisation time (cf. eq. [7|). Therefore, if 



>fv 



(15) 



then the giant embryo may collapse as a whole. Detailed 
calcu lations of this proce ss in the planet formation context 
(e.g.. lBodenheimerlll974l l show that a Jupiter mass embryo 
only collapses in this way after ~ 10'' years, which appears 
longer than the migration time for "reasonable" parameters. 
However, more massive embryos start off hotter and should 
reach the second collapse stage sooner (see paper II). Thus 
it is possible that giant planets more massive than Jupiter 
may have no solid cores. 



5 ASTROPHYSICAL IMPLICATIONS 

When does planet formation stop? The view advanced 
in ij tl4. 1.1 1 and 14.1.31 is an extreme version of the tidal down- 
sizing hypothesis, where everything "interesting" occurs in- 
side the GEs before they are disrupted. This scenario is un- 
likely to be realistic. When the GEs are dispersed, their rem- 
nants are still embedded in a protoplanetary disc. The disc 
may be less massive, but as observations show it may persist 
for a few million years. Therefore, it is likely that the planet 
building action continues at that stage, now as in the CA 
model - solid cores may grow further by accreting smaller 
solid debris; giant planets may accrete more gas and also mi- 
grate radially. Therefore, the most reasonable point of view 
is that the tidal downsizing hypothesis, if correct, may be a 
new way to begin planetar y system formation process (see 
also IClarke fc Lodatol |2009| ) . 

Metalicity correlation. Observations show that the frac- 
tion of stars orbited by a giant planet increases s trongly 
with metalicity of the star (jFischer fc Valentill2005^ . As in 
the CA model, more solids in the gas enable more massive 
solid cores to be built. We also assume that if no solid core 
is built, the GE is completely destroyed close to the star; all 
of the GE material presumably accretes onto the star. 

In detail, core formation requires igr < tvap (32IU- Dust 
opacity, k, , is likely to be proportional to the metalicity 
and the grain mass fraction, fg. In such a model, the grain 
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growth timescale behaves as igr oc /, 



-5/21 



(eq. [HI), whereas 



ivap oc fg. A higher fg thus gives the GEmore time to build 
a solid core. Further, iMeru fc Batd (|2010l ) find that lower 
dust opacities discs fragment closer to the star, reducing 
migration time (cf. eq. [5]). Low metalicity GEs are likely to 
vaporise their grains before they make a solid core, and they 
also migrate inward too quickly to their ultimate demise. 

Challenges and open questions. "Tidal downsizing" is 
a complex hypothesis that may turn out to not work for at 
least the following reasons. (1) giant planet embryos need to 
be isolated or slowly accreting. If they pile up mass quickly 
by accretion from the disc, they become too hot to support 
grain sedimentation (paper I). (2) We assumed that GEs mi- 
grate as planets do. This is probably a good approximation 
when the radius of the embryo -Rcmb ^ -Rh, but much more 
work is needed to confirm this. (3) We have seen that the 
timescales of the important processes (see H4.ip do depend 
strongly on dust opacity and other parameters. Real proto- 
planetary discs may be such that grain growth inside the 
GEs is precluded or the embryos are disrupted too quickly. 



6 CONCLUSIONS 

The tidal downsizing hypothesis combines some of the planet 
building processes from both the CA and the GI models, pre- 
serving the strengths but not the weaknesses of these theo- 
ries. Open questions remain, but on the balance we feel that 
the hypothesis deserves a further detailed consideration by 
the planet formation community. We note that, having sub- 
mitted this paper and paper I, we discovered an independent 
suggestion on the possibility of terr estrial core format ion by 
tidally disrupted giant embryos bv lBolev et al.l ((20101). 
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APPENDIX: IRRADIATIVE EVAPORATION 

Irradiati on of the GE by th e star can affect the former 
strongly (|Cameron et al.lll982h . First we estimate the radius 
i^iock below which the temperature of the surface layers of 
the GE is set by the irradiation rather than the internal 
radiation. This is found by equating the stellar radiation lu- 
minosity incident on the GE, Lirr = I/*7r_R^^ij/(47r_Rp), to 
the isolated GE radiative luminosity, Lf^{t). Since Lfc{t) oc 
T{t)-^ oc (l + 2f/fo)"^''^ for the GE (paper I), the result for 
-Riock at i = tii and L, — Lq is 



r> ir, ATT -l;3/4 l/4„-3/8 

-Riock ~ 12 AU nil kj a_-^ itj 



H 



0.2R^ 



1/2 



Irradiation can also disrupt the GE when the irradiation 
temperature. Tin- — (Lir,/47r_R^jnj,o-)^''* is comparable to the 
mean temperature (equation [1} inside the GE, e.g., Tlrr = 
rievT (t), with r]sv ^ 1 . For one of the fixed planet location 
runs I Cameron et al.l l|l982h found -qev ~ 0.5 (cf. their Fig. 
2). As Rev = {L,/'iTvaT^y^^{to/2tii), we have 



Re 



1 AU k,a-iK. 



-3/2 



H 



0.2i?, 



0.5 



-1/3 



(17) 



(16) 



For fiducial parameters, tidal disruption of GEs occurs ear- 
lier (at larger radii) than the irradiative evaporation. Fur- 
thermore, as noted by the referee, the giant embryo may be 
shadowed by the disc if the disc scale-height, H, is larger 
than the GE size, Rcmh- Thus we expect tidal disruption of 
the GEs to be more frequent. Note however different param- 
eter dependencies in ea. ll2l and ll7l In general both processes 
may be important in removing GE envelopes, with irradia- 
tion becoming more important for higher mass stars. 
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